Computational study of diffraction image formation from XFEL irradiated single ribosome molecule.

Single particle imaging at atomic resolution is perhaps one of the most desired goals for ultrafast X-ray science with X-ray free-electron lasers. Such a capability would create great opportunity within the biological sciences, as high-resolution structural information of biosamples that may not crystallize is essential for many research areas therein. In this paper, we report on a comprehensive computational study of diffraction image formation during single particle imaging of a macromolecule, containing over one hundred thousand non-hydrogen atoms. For this study, we use a dedicated simulation framework, SIMEX, available at the European XFEL facility. Our results demonstrate the full feasibility of computational single-particle imaging studies for biological samples of realistic size. This finding is important as it shows that the SIMEX platform can be used for simulations to inform relevant single-particle-imaging experiments and help to establish optimal parameters for these experiments. This will enable more focused and more efficient single-particle-imaging experiments at XFEL facilities, making the best use of the resource-intensive XFEL operation.

SARS-CoV-2 Mpro responds to oxidation by forming disulfide and NOS/SONOS bonds.

The main protease (Mpro) of SARS-CoV-2 is critical for viral function and a key drug target. Mpro is only active when reduced; turnover ceases upon oxidation but is restored by re-reduction. This suggests the system has evolved to survive periods in an oxidative environment, but the mechanism of this protection has not been confirmed. Here, we report a crystal structure of oxidized Mpro showing a disulfide bond between the active site cysteine, C145, and a distal cysteine, C117. Previous work proposed this disulfide provides the mechanism of protection from irreversible oxidation. Mpro forms an obligate homodimer, and the C117-C145 structure shows disruption of interactions bridging the dimer interface, implying a correlation between oxidation and dimerization. We confirm dimer stability is weakened in solution upon oxidation. Finally, we observe the protein's crystallization behavior is linked to its redox state. Oxidized Mpro spontaneously forms a distinct, more loosely packed lattice. Seeding with crystals of this lattice yields a structure with an oxidation pattern incorporating one cysteine-lysine-cysteine (SONOS) and two lysine-cysteine (NOS) bridges. These structures further our understanding of the oxidative regulation of Mpro and the crystallization conditions necessary to study this structurally.

MolDStruct: Modeling the dynamics and structure of matter exposed to ultrafast x-ray lasers with hybrid collisional-radiative/molecular dynamics.

We describe a method to compute photon-matter interaction and atomic dynamics with x-ray lasers using a hybrid code based on classical molecular dynamics and collisional-radiative calculations. The forces between the atoms are dynamically determined based on changes to their electronic occupations and the formation of a free electron cloud created from the irradiation of photons in the x-ray spectrum. The rapid transition from neutral solid matter to dense plasma phase allows the use of screened potentials, reducing the number of non-bonded interactions. In combination with parallelization through domain decomposition, the hybrid code handles large-scale molecular dynamics and ionization. This method is applicable for large enough samples (solids, liquids, proteins, viruses, atomic clusters, and crystals) that, when exposed to an x-ray laser pulse, turn into a plasma in the first few femtoseconds of the interaction. We present four examples demonstrating the applicability of the method. We investigate the non-thermal heating and scattering of bulk water and damage-induced dynamics of a protein crystal using an x-ray pump-probe scheme. In both cases, we compare to the experimental data. For single particle imaging, we simulate the ultrafast dynamics of a methane cluster exposed to a femtosecond x-ray laser. In the context of coherent diffractive imaging, we study the fragmentation as given by an x-ray pump-probe setup to understand the evolution of radiation damage in the time range of hundreds of femtoseconds.

Observation of a single protein by ultrafast X-ray diffraction.

The idea of using ultrashort X-ray pulses to obtain images of single proteins frozen in time has fascinated and inspired many. It was one of the arguments for building X-ray free-electron lasers. According to theory, the extremely intense pulses provide sufficient signal to dispense with using crystals as an amplifier, and the ultrashort pulse duration permits capturing the diffraction data before the sample inevitably explodes. This was first demonstrated on biological samples a decade ago on the giant mimivirus. Since then, a large collaboration has been pushing the limit of the smallest sample that can be imaged. The ability to capture snapshots on the timescale of atomic vibrations, while keeping the sample at room temperature, may allow probing the entire conformational phase space of macromolecules. Here we show the first observation of an X-ray diffraction pattern from a single protein, that of Escherichia coli GroEL which at 14 nm in diameter is the smallest biological sample ever imaged by X-rays, and demonstrate that the concept of diffraction before destruction extends to single proteins. From the pattern, it is possible to determine the approximate orientation of the protein. Our experiment demonstrates the feasibility of ultrafast imaging of single proteins, opening the way to single-molecule time-resolved studies on the femtosecond timescale.

Revealing the origins of vortex cavitation in a Venturi tube by high speed X-ray imaging.

Hydrodynamic cavitation is useful in many processing applications, for example, in chemical reactors, water treatment and biochemical engineering. An important type of hydrodynamic cavitation that occurs in a Venturi tube is vortex cavitation known to cause luminescence whose intensity is closely related to the size and number of cavitation events. However, the mechanistic origins of bubbles constituting vortex cavitation remains unclear, although it has been concluded that the pressure fields generated by the cavitation collapse strongly depends on the bubble geometry. The common view is that vortex cavitation consists of numerous small spherical bubbles. In the present paper, aspects of vortex cavitation arising in a Venturi tube were visualized using high-speed X-ray imaging at SPring-8 and European XFEL. It was discovered that vortex cavitation in a Venturi tube consisted of angulated rather than spherical bubbles. The tangential velocity of the surface of vortex cavitation was assessed considering the Rankine vortex model.

Structure of the Lysinibacillus sphaericus Tpp49Aa1 pesticidal protein elucidated from natural crystals using MHz-SFX.

The Lysinibacillus sphaericus proteins Tpp49Aa1 and Cry48Aa1 can together act as a toxin toward the mosquito Culex quinquefasciatus and have potential use in biocontrol. Given that proteins with sequence homology to the individual proteins can have activity alone against other insect species, the structure of Tpp49Aa1 was solved in order to understand this protein more fully and inform the design of improved biopesticides. Tpp49Aa1 is naturally expressed as a crystalline inclusion within the host bacterium, and MHz serial femtosecond crystallography using the novel nanofocus option at an X-ray free electron laser allowed rapid and high-quality data collection to determine the structure of Tpp49Aa1 at 1.62 Å resolution. This revealed the packing of Tpp49Aa1 within these natural nanocrystals as a homodimer with a large intermolecular interface. Complementary experiments conducted at varied pH also enabled investigation of the early structural events leading up to the dissolution of natural Tpp49Aa1 crystals-a crucial step in its mechanism of action. To better understand the cooperation between the two proteins, assays were performed on a range of different mosquito cell lines using both individual proteins and mixtures of the two. Finally, bioassays demonstrated Tpp49Aa1/Cry48Aa1 susceptibility of Anopheles stephensi, Aedes albopictus, and Culex tarsalis larvae-substantially increasing the potential use of this binary toxin in mosquito control.

A phenomenological model of the X-ray pulse statistics of a high-repetition-rate X-ray free-electron laser.

Many coherent imaging applications that utilize ultrafast X-ray free-electron laser (XFEL) radiation pulses are highly sensitive to fluctuations in the shot-to-shot statistical properties of the source. Understanding and modelling these fluctuations are key to successful experiment planning and necessary to maximize the potential of XFEL facilities. Current models of XFEL radiation and their shot-to-shot statistics are based on theoretical descriptions of the source and are limited in their ability to capture the shot-to-shot intensity fluctuations observed experimentally. The lack of accurate temporal statistics in simulations that utilize these models is a significant barrier to optimizing and interpreting data from XFEL coherent diffraction experiments. Presented here is a phenomenological model of XFEL radiation that is capable of capturing the shot-to-shot statistics observed experimentally using a simple time-dependent approximation of the pulse wavefront. The model is applied to reproduce non-stationary shot-to-shot intensity fluctuations observed at the European XFEL, whilst accurately representing the single-shot properties predicted by FEL theory. Compared with previous models, this approach provides a simple, robust and computationally inexpensive method of generating statistical representations of XFEL radiation.

Form factor determination of biological molecules with X-ray free electron laser small-angle scattering (XFEL-SAS).

Free-electron lasers (FEL) are revolutionizing X-ray-based structural biology methods. While protein crystallography is already routinely performed at FELs, Small Angle X-ray Scattering (SAXS) studies of biological macromolecules are not as prevalent. SAXS allows the study of the shape and overall structure of proteins and nucleic acids in solution, in a quasi-native environment. In solution, chemical and biophysical parameters that have an influence on the structure and dynamics of molecules can be varied and their effect on conformational changes can be monitored in time-resolved XFEL and SAXS experiments. We report here the collection of scattering form factors of proteins in solution using FEL X-rays. The form factors correspond to the scattering signal of the protein ensemble alone; the scattering contributions from the solvent and the instrument are separately measured and accurately subtracted. The experiment was done using a liquid jet for sample delivery. These results pave the way for time-resolved studies and measurements from dilute samples, capitalizing on the intense and short FEL X-ray pulses.

Water layer and radiation damage effects on the orientation recovery of proteins in single-particle imaging at an X-ray free-electron laser.

The noise caused by sample heterogeneity (including sample solvent) has been identified as one of the determinant factors for a successful X-ray single-particle imaging experiment. It influences both the radiation damage process that occurs during illumination as well as the scattering patterns captured by the detector. Here, we investigate the impact of water layer thickness and radiation damage on orientation recovery from diffraction patterns of the nitrogenase iron protein. Orientation recovery is a critical step for single-particle imaging. It enables to sort a set of diffraction patterns scattered by identical particles placed at unknown orientations and assemble them into a 3D reciprocal space volume. The recovery quality is characterized by a "disconcurrence" metric. Our results show that while a water layer mitigates protein damage, the noise generated by the scattering from it can introduce challenges for orientation recovery and is anticipated to cause problems in the phase retrieval process to extract the desired protein structure. Compared to these disadvantageous effects due to the thick water layer, the effects of radiation damage on the orientation recovery are relatively small. Therefore, minimizing the amount of residual sample solvent should be considered a crucial step in improving the fidelity and resolution of X-ray single-particle imaging experiments.

Online dynamic flat-field correction for MHz microscopy data at European XFEL.

The high pulse intensity and repetition rate of the European X-ray Free-Electron Laser (EuXFEL) provide superior temporal resolution compared with other X-ray sources. In combination with MHz X-ray microscopy techniques, it offers a unique opportunity to achieve superior contrast and spatial resolution in applications demanding high temporal resolution. In both live visualization and offline data analysis for microscopy experiments, baseline normalization is essential for further processing steps such as phase retrieval and modal decomposition. In addition, access to normalized projections during data acquisition can play an important role in decision-making and improve the quality of the data. However, the stochastic nature of X-ray free-electron laser sources hinders the use of standard flat-field normalization methods during MHz X-ray microscopy experiments. Here, an online (i.e. near real-time) dynamic flat-field correction method based on principal component analysis of dynamically evolving flat-field images is presented. The method is used for the normalization of individual X-ray projections and has been implemented as a near real-time analysis tool at the Single Particles, Clusters, and Biomolecules and Serial Femtosecond Crystallography (SPB/SFX) instrument of EuXFEL.

Heterologous expression, purification and structural features of native Dictyostelium discoideum dye-decolorizing peroxidase bound to a natively incorporated heme.

The Dictyostelium discoideum dye-decolorizing peroxidase (DdDyP) is a newly discovered peroxidase, which belongs to a unique class of heme peroxidase family that lacks homology to the known members of plant peroxidase superfamily. DdDyP catalyzes the H2O2-dependent oxidation of a wide-spectrum of substrates ranging from polycyclic dyes to lignin biomass, holding promise for potential industrial and biotechnological applications. To study the molecular mechanism of DdDyP, highly pure and functional protein with a natively incorporated heme is required, however, obtaining a functional DyP-type peroxidase with a natively bound heme is challenging and often requires addition of expensive biosynthesis precursors. Alternatively, a heme in vitro reconstitution approach followed by a chromatographic purification step to remove the excess heme is often used. Here, we show that expressing the DdDyP peroxidase in ×2 YT enriched medium at low temperature (20°C), without adding heme supplement or biosynthetic precursors, allows for a correct native incorporation of heme into the apo-protein, giving rise to a stable protein with a strong Soret peak at 402 nm. Further, we crystallized and determined the native structure of DdDyP at a resolution of 1.95 Å, which verifies the correct heme binding and its geometry. The structural analysis also reveals a binding of two water molecules at the distal site of heme plane bridging the catalytic residues (Arg239 and Asp149) of the GXXDG motif to the heme-Fe(III) via hydrogen bonds. Our results provide new insights into the geometry of native DdDyP active site and its implication on DyP catalysis.

XFEL Microcrystallography of Self-Assembling Silver n-Alkanethiolates.

New synthetic hybrid materials and their increasing complexity have placed growing demands on crystal growth for single-crystal X-ray diffraction analysis. Unfortunately, not all chemical systems are conducive to the isolation of single crystals for traditional characterization. Here, small-molecule serial femtosecond crystallography (smSFX) at atomic resolution (0.833 Å) is employed to characterize microcrystalline silver n-alkanethiolates with various alkyl chain lengths at X-ray free electron laser facilities, resolving long-standing controversies regarding the atomic connectivity and odd-even effects of layer stacking. smSFX provides high-quality crystal structures directly from the powder of the true unknowns, a capability that is particularly useful for systems having notoriously small or defective crystals. We present crystal structures of silver n-butanethiolate (C4), silver n-hexanethiolate (C6), and silver n-nonanethiolate (C9). We show that an odd-even effect originates from the orientation of the terminal methyl group and its role in packing efficiency. We also propose a secondary odd-even effect involving multiple mosaic blocks in the crystals containing even-numbered chains, identified by selected-area electron diffraction measurements. We conclude with a discussion of the merits of the synthetic preparation for the preparation of microdiffraction specimens and compare the long-range order in these crystals to that of self-assembled monolayers.

Hard X-ray stereographic microscopy for single-shot differential phase imaging.

The characterisation of fast phenomena at the microscopic scale is required for the understanding of catastrophic responses of materials to loads and shocks, the processing of materials by optical or mechanical means, the processes involved in many key technologies such as additive manufacturing and microfluidics, and the mixing of fuels in combustion. Such processes are usually stochastic in nature and occur within the opaque interior volumes of materials or samples, with complex dynamics that evolve in all three dimensions at speeds exceeding many meters per second. There is therefore a need for the ability to record three-dimensional X-ray movies of irreversible processes with resolutions of micrometers and frame rates of microseconds. Here we demonstrate a method to achieve this by recording a stereo phase-contrast image pair in a single exposure. The two images are combined computationally to reconstruct a 3D model of the object. The method is extendable to more than two simultaneous views. When combined with megahertz pulse trains of X-ray free-electron lasers (XFELs) it will be possible to create movies able to resolve 3D trajectories with velocities of kilometers per second.

Automatic bad-pixel mask maker for X-ray pixel detectors with application to serial crystallography.

X-ray crystallography has witnessed a massive development over the past decade, driven by large increases in the intensity and brightness of X-ray sources and enabled by employing high-frame-rate X-ray detectors. The analysis of large data sets is done via automatic algorithms that are vulnerable to imperfections in the detector and noise inherent with the detection process. By improving the model of the behaviour of the detector, data can be analysed more reliably and data storage costs can be significantly reduced. One major requirement is a software mask that identifies defective pixels in diffraction frames. This paper introduces a methodology and program based upon concepts of machine learning, called robust mask maker (RMM), for the generation of bad-pixel masks for large-area X-ray pixel detectors based on modern robust statistics. It is proposed to discriminate normally behaving pixels from abnormal pixels by analysing routine measurements made with and without X-ray illumination. Analysis software typically uses a Bragg peak finder to detect Bragg peaks and an indexing method to detect crystal lattices among those peaks. Without proper masking of the bad pixels, peak finding methods often confuse the abnormal values of bad pixels in a pattern with true Bragg peaks and flag such patterns as useful regardless, leading to storage of enormous uninformative data sets. Also, it is computationally very expensive for indexing methods to search for crystal lattices among false peaks and the solution may be biased. This paper shows how RMM vastly improves peak finders and prevents them from labelling bad pixels as Bragg peaks, by demonstrating its effectiveness on several serial crystallography data sets.

Pump-probe capabilities at the SPB/SFX instrument of the European XFEL.

Pump-probe experiments at X-ray free-electron laser (XFEL) facilities are a powerful tool for studying dynamics at ultrafast and longer timescales. Observing the dynamics in diverse scientific cases requires optical laser systems with a wide range of wavelength, flexible pulse sequences and different pulse durations, especially in the pump source. Here, the pump-probe instrumentation available for measurements at the Single Particles, Clusters, and Biomolecules and Serial Femtosecond Crystallography (SPB/SFX) instrument of the European XFEL is reported. The temporal and spatial stability of this instrumentation is also presented.

Towards real-time analysis of liquid jet alignment in serial femtosecond crystallography.

Liquid sample delivery systems are used extensively for serial femtosecond crystallography at X-ray free-electron lasers (XFELs). However, misalignment of the liquid jet and the XFEL beam leads to the X-rays either partially or completely missing the sample, resulting in sample wastage and a loss of experiment time. Implemented here is an algorithm to analyse optical images using machine vision to determine whether there is overlap of the X-ray beam and liquid jet. The long-term goal is to use the output from this algorithm to implement an automated feedback mechanism to maintain constant alignment of the X-ray beam and liquid jet. The key elements of this jet alignment algorithm are discussed and its performance is characterized by comparing the results with a manual analysis of the optical image data. The success rate of the algorithm for correctly identifying hits is quantified via a similarity metric, the Dice coefficient. In total four different nozzle designs were used in this study, yielding an overall Dice coefficient of 0.98.

Shot-to-shot two-dimensional photon intensity diagnostics within megahertz pulse-trains at the European XFEL.

Characterizing the properties of X-ray free-electron laser (XFEL) sources is a critical step for optimization of performance and experiment planning. The recent availability of MHz XFELs has opened up a range of new opportunities for novel experiments but also highlighted the need for systematic measurements of the source properties. Here, MHz-enabled beam imaging diagnostics developed for the SPB/SFX instrument at the European XFEL are exploited to measure the shot-to-shot intensity statistics of X-ray pulses. The ability to record pulse-integrated two-dimensional transverse intensity measurements at multiple planes along an XFEL beamline at MHz rates yields an improved understanding of the shot-to-shot photon beam intensity variations. These variations can play a critical role, for example, in determining the outcome of single-particle imaging experiments and other experiments that are sensitive to the transverse profile of the incident beam. It is observed that shot-to-shot variations in the statistical properties of a recorded ensemble of radiant intensity distributions are sensitive to changes in electron beam current density. These changes typically occur during pulse-distribution to the instrument and are currently not accounted for by the existing suite of imaging diagnostics. Modulations of the electron beam orbit in the accelerator are observed to induce a time-dependence in the statistics of individual pulses - this is demonstrated by applying radio-frequency trajectory tilts to electron bunch-trains delivered to the instrument. We discuss how these modifications of the beam trajectory might be used to modify the statistical properties of the source and potential future applications.

Tree-Code Based Improvement of Computational Performance of the X-ray-Matter-Interaction Simulation Tool XMDYN.

In this work, we report on incorporating for the first time tree-algorithm based solvers into the molecular dynamics code, XMDYN. XMDYN was developed to describe the interaction of ultrafast X-ray pulses with atomic assemblies. It is also a part of the simulation platform, SIMEX, developed for computational single-particle imaging studies at the SPB/SFX instrument of the European XFEL facility. In order to improve the XMDYN performance, we incorporated the existing tree-algorithm based Coulomb solver, PEPC, into the code, and developed a dedicated tree-algorithm based secondary ionization solver, now also included in the XMDYN code. These extensions enable computationally efficient simulations of X-ray irradiated large atomic assemblies, e.g., large protein systems or viruses that are of strong interest for ultrafast X-ray science. The XMDYN-based preparatory simulations can now guide future single-particle-imaging experiments at the free-electron-laser facility, EuXFEL.

A multi-million image Serial Femtosecond Crystallography dataset collected at the European XFEL.

Serial femtosecond crystallography is a rapidly developing method for determining the structure of biomolecules for samples which have proven challenging with conventional X-ray crystallography, such as for membrane proteins and microcrystals, or for time-resolved studies. The European XFEL, the first high repetition rate hard X-ray free electron laser, provides the ability to record diffraction data at more than an order of magnitude faster than previously achievable, putting increased demand on sample delivery and data processing. This work describes a publicly available serial femtosecond crystallography dataset collected at the SPB/SFX instrument at the European XFEL. This dataset contains information suitable for algorithmic development for detector calibration, image classification and structure determination, as well as testing and training for future users of the European XFEL and other XFELs.

Unsupervised learning approaches to characterizing heterogeneous samples using X-ray single-particle imaging.

One of the outstanding analytical problems in X-ray single-particle imaging (SPI) is the classification of structural heterogeneity, which is especially difficult given the low signal-to-noise ratios of individual patterns and the fact that even identical objects can yield patterns that vary greatly when orientation is taken into consideration. Proposed here are two methods which explicitly account for this orientation-induced variation and can robustly determine the structural landscape of a sample ensemble. The first, termed common-line principal component analysis (PCA), provides a rough classification which is essentially parameter free and can be run automatically on any SPI dataset. The second method, utilizing variation auto-encoders (VAEs), can generate 3D structures of the objects at any point in the structural landscape. Both these methods are implemented in combination with the noise-tolerant expand-maximize-compress (EMC) algorithm and its utility is demonstrated by applying it to an experimental dataset from gold nanoparticles with only a few thousand photons per pattern. Both discrete structural classes and continuous deformations are recovered. These developments diverge from previous approaches of extracting reproducible subsets of patterns from a dataset and open up the possibility of moving beyond the study of homogeneous sample sets to addressing open questions on topics such as nanocrystal growth and dynamics, as well as phase transitions which have not been externally triggered.